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Abstract 

We describe an algorithm for testing the completeness of caps in 
PG(r, q), q even. It allowed us to check that the 95256-cap in PG(12,4) 
recently found by Fu el al. (see m) is complete. 


1 Introduction 


Let PG(r, q) be the r-dimensional projective space over the Galois field ¥q. 
An n-cap in PG(r, q) is a set of points no three of which are collinear. An 
n-cap in PG(r, g) is called complete if it is not contained in an (n + l)-cap 
in PG(r, g); see [T6] . 

The points of a complete n-cap in PG(r—1, g) can be treated as columns of 
a parity check matrix of an [n,n — r, 4]g linear code with the exceptions of the 
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complete 5-cap in PG(3, 2) and the complete 11-cap in PG(4, 3) correspond¬ 
ing to the binary [5,1, 5 ]2 code and to the Golay [11, 6 , 5 ]3 code respectively. 

An n-cap in PG(r, q) of maximal size is called a maximal cap in PG(r, q). 
A classical problem on caps is to determine the maximal size of complete caps 
in PG(r, q). This is also known as the packing problem; see [1^. Denote the 
size of a maximal cap in PG(r, q) as m 2 {r, q), and the largest size of a known 
complete cap as m 2 (r, q). 

Of particular interest is the case g = 4, due the connection with quan¬ 
tum error correction established in [9], where a class of quantum codes, the 
quantum stabilizer codes, is described in terms of certain additive quaternary 
codes. 

Additive quaternary codes are dehned over F 4 but are linear over F 2 . If 
we restrict considering quaternary quantum codes that are indeed F 4 -hnear 
then we have the following definition; see IHEl: 

Definition 1.1. A linear quaternary quantum stabilizer code is a subspace 
C C F4 such that C C where duality is with respect to the Hermitian 

inner product. 

Here the Hermitian inner product of a; = (xi, ...,Xn) and y = {yi, ...,yn) 
is {x,y) = where y = y"^. The reason for this definition is that 

a linear quaternary quantum stabilizer code C of length tt,, dimension r and 
dual distance > d (equivalently: of strength > d) allows the construction of 
a pure quantum stabilizer code [[n,n — 2r,c[\]i] see [3 Theorem 1]. 

A pure quantum code [[n,n — 2r, 4]] which is linear over F 4 is obtained 
from a cap satisfying certain conditions; see [H Theorem 2.8]: 

Definition 1.2. A cap ^ in PG{r — 1,4) is a quantum cap if it is not con¬ 
tained in a hyperplane and if it satisfies the following equivalent conditions: 

• each hyperplane meets the cap in the same parity as the cardinality of 
the cap; 

• the corresponding quaternary [n,r] 4 ^-code has all its weights even; 

• the corresponding quaternary [n, r] 4 -code is self-orthogonal with respect 
to the Hermitian inner product. 

Theorem 1.3. The following are equivalent: 

• A pure stabilizer quantum code [[n, n — 2r, 4]] which is linear over ¥ 4 ^. 

• A quantum n-cap in PG{r — 1,4). 
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Much work on caps has been done, see [3-7,10-21], The value of m 2 (r, 4) 
is known for k < A\ 7712 ( 2 ,4) = 6, ?t 72(3,4) = 17, and 7772(4,4) = 41. 

In [H] it is proved that 7772 ( 8 , 4 ) = 2136, 7772 ( 9 ,4) = 5124, 7772(10,4) = 
15840, 7772 ( 11 , 4 ) = 36084 and they also give a 95256-cap in PG(12,4). 

Their results have been obtained by computer-supported recursive con¬ 
structions. They also present an algorithm for checking completeness of a 
cap based on a bijective map between points in PG{r, 4) and a subset / of 
the positive integer set N. 

This algorithm allowed checking the completeness of the caps for k < 
11, but it is too computationally expensive for the case k = 12. As they 
wrote: “But as for checking completeness of larger caps in PG(r, 4), r > 12, 
new algorithms are needed.”; see Cl Section 5]. We propose a new fast 
algorithm that allowed to face also this case: we verihed that the 95256-cap 
in PG(12,4) is complete, so 7772 ( 12 , 4 ) = 95256. Our algorithm is based on a 
compact representation of the points of PG(r, q), q even, and on minimizing 
the computational costs of the operations more often performed during the 
check of the completeness of the cap. 

Section 2 describes the algorithm and applies it in PG(12,4). Section 3 
contains the generalization of the algorithm to other even values of q and 
other dimensions. 


2 A new algorithm for checking completeness 
of a cap 

In [H] an algorithm for checking completeness of a cap ^ in PG{r, 4) is 
presented. It is based on a bijective map 0 between points in PG{r, 4) and 
a subset T(r) of the positive integer set N: 

0 : PG(r,4) ^ T(r), 

0:P^0(P), 


where 


P = {Xo,Xi, ...,Xr)'^, 

(j){P) = P^xo + 4’-ixi+-+4p-i + P. 
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Table 1; The time and space cost of the algorithm of [13] 


Size of cap 

2136 

5124 

15840 

36084 

Time 

7805 

46244 

428029 

2261301 


It can be easily seen that a cap is complete if and only if each point P 
of PG(r, q) not belonging to the cap lies on a secant line of the cap. In this 
case we say that P is covered. 

To keep track of the covering of the points, a vector U of size |T(r)| is 
used. Initially all elements of U are set to be 1. 

Then all pairs of points of the cap are considered. For each pair of points 
(Pj, Pj) the three other points belonging to the line through Pj, Pj are com¬ 
puted. To do this all linear combination Q = aPi -|- P 2 , a G F 4 \{ 0 } are 
computed. The point Q is normalized, choosing a representation with the 
leftmost non-zero coordinate equal to 1. Finally the position 0((5) of U is 
set to 0. The process continues until all elements of U became 0 or all pairs 
of points of the cap have been considered. 

At the end the cap is complete if and only if all elements of U are 0. 
Table 1 reports the time cost of the algorithm using an Intel(R) Xeon(R) 
CPU E5504 @ 2.00GHz; see [HI Table 1]. However, the paper does not 
mention the unit of time used in Table 1. 

We devised a new algorithm for checking the completeness of a cap in 
PG{r, 4) choosing a representation that optimizes the computational cost 
of the main operations of the previous algorithm: the computation of Q = 
aPi + P 2 and the normalization of a point Q. 

Let be F 4 = {0, 1 , 0 ;,^}, where = U, U = u + 1, and = 1. If we 
dehne a representation function p : F 4 ^ N as in the following: 

p(0) —)■ 0, p(l) —)■ 1, p{u) = 2, p{uJ) = 3, 


then we have 

a + b = p{a)''p{b), a, fe G F4, 

where " is the bitwise exclusive or operator. 

Moreover, if P = {xo,Xi, ...,Xr)'^ then the binary representation of 0(P) 
is p{xo)p{xi). ..p{Xr). 

It means that if Q = Pi -|- P 2 then (j){Q) = 0(Pi)'0(P2). 
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This allows the computation of the sum of two points of PG{r, 4) by one 
integer operation. 

The multiplication of one point P by a scalar is applied only to the 
points of the cap. It can be pre-computed before the beginning of the check 
for completeness, so at the cost of having a data structure of size 3|^| all 
multiplications by a scalar are avoided. 

The other expensive operation of the algorithm in [IT] is the normalization 
of a point P. It should be to computed each time Q = «Pi + P 2 , Pi, P 2 G is 
computed, i.e. 3/2|^p times. We propose a trade-off between computational 
time and memory space: we use a vector U of size 3|PG(r, 4)| to keep trace 
of the fact that a point Q is covered by or not; initially all elements of U 
are set equal to 0. When Q = aPi + P 2 is computed, then the element 0(Q) 
is set equal to 1 without before normalizing Q. In this way all the 3/2|‘^p 
normalization operations are avoided. At the end, when the covering of all 
points of PG(r, 4) is tested, first the normalized form of a point Q is tested 
checking the element of U of position if it is not covered also (j){u}Q) 

and 4>{lJQ) are checked: if any of 0(Q), 0(a;Q), 0(aJQ) is equal to 1, then Q 
is covered. 

Let be n = 1^1, m = |PG(r, q)\, i = the size of an integer, c = the size of 
a character. The total cost of our algorithm is: 

space: 3n ■ i + 3m ■ c + Ci, 
time: 3c2n -|- 3/2c2n^ -l- dcsm ■ c + c^, 

where ci,..., C 4 are constants. 

The algorithm has been implemented in C language. 

Table 2 reports the time and space cost of the algorithm using an Intel(R) 
Core(TM) i7-4510U CPU @ 2.00GHz; space is measured in megabyte, while 
time is measured in milliseconds. 

We tested the completeness of the caps presented in [Tl] in PG{r,q), 
r = 8 ,..., 12. Note that for constructing the 5124-cap in PG(9,4) and 
the 36084-cap in PG(11,4) we were not able to obtain a cap following the 
selection of columns suggested in [TTl Section 3.2]. For the 5124-cap we had 
to exclude the vectors for j e {2,14, 24, 25} instead of j G {2,14,15, 24} 
as suggested in the paper, whereas for the 36084-cap we had to exclude the 
vectors for i G {1,..., 16, 271} instead of i G {1, 256,..., 271} as suggested 
in the paper. 
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Table 2: The time and space cost of the new algorithm 


Size of cap 

2136 

5124 

15840 

36084 

95256 

Space (Mb) 

5 

18 

76 

324 

1382 

Time (Milliseconds) 

32 

78 

707 

3574 

98321 


3 Generalization of the algorithm 

In the previous section we applied our algorithm in PG(12,4). The key idea, 
choosing a representation for the elements of and the points of PG(r, q) 
that minimize the computational cost of the operations most often performed 
during the test of completeness of a cap, can be applied for every even q. 

When testing the completeness of a cap ^ the value Q = aPi + P 2 , 
Pi, P 2 & ^Oi E Fg\{0} has to be computed. 

To avoid to compute the same value aP several times, it is convenient to 
compute it at the beginning of the algorithm and store the results. There¬ 
fore the main operation to compute remains the sum between two vectors 
representing points of PG(r, q). 

When considering a representation of F 2 fc, we can either choose a form 
that facilitate the computation of multiplication (we see the non-zero ele¬ 
ments of F 2 fe as powers of the primitive element) or can choose a form that 
facilitate the computation of addition (we see the the elements of Fg as poly¬ 
nomials of F 2 [X] of degree less than k] addition is dehned in the natural way, 
whereas multiplication is dehned modulo a hxed irreducible polynomial of 
degree k). 

We choose the latter representation and dehne p : F 2 fc —)■ {0,..., 2^ — 1} 
as p : p{x) I—)■ p(2). We have that p{p{x) + q{x)) = p{p{x))''p{q{x)), where 
is the bitwise exclusive or. It means that in this representation addition on 
F 2 fc reduces to one bitwise arithmetic operation on integers. 

Moreover we can dehne a representation of the points of PG{r, 2^) in the 
following way: 


0 : PG{r,2^) N, 
0:P^0(P), 


where 
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P = {Xo,Xi, ...,Xr)'^, 

0(P) = (2^)V(xo) + (2^)'-V(xi) + • • • + (2^)p(p-i) + p(P). 

In this way a point P of PG{r, 2^) is represented by an integer n. If 
we consider the binary representation of n, the coordinate Xj is represented 
by the bits of n in position {r — i)k + 1.. .{r — i + l)k that are the binary 
representation p(xj). To compnte 0(Q), Q = P 1 +P 2 h is snfficient compnting 
0(Pi)''0(P2), just one bitwise arithmetic operation. In a real implementation, 
usually an (unsigned) integer has a 32 bit representation, so 0(P) can be 
represented by a single integer if kr < 32, otherwise more integers are needed. 

Our algorithm trades computational time for memory space. Let be n = 
I'^l, m = |PG(r, 2^)1, we need ?7,(2^ —1) integers to represent aP, P E a & 
F2fc\{0} and m(2^ —1) booleans to represent the fact that aP, P G PG{r, 2^), 
a G F2fe\{0} is saturated or not. As n < m, the latter value is more relevant. 
If the memory requested by the algorithm is too big, then memory space can 
be traded for computational time. For example the set of points of PG{r, 2^) 
can be divided into s subsets small enough to be represented. Than the test 
for completeness can be repeated s times; each time the completeness of the 
points of one subset is tested. Note that the computations for the different 
subsets are independent, so they can be performed in parallel. This is a form 
of parallelism based on the splitting of data: it is simple and effective. 

4 Conclusion 

We presented an algorithm for checking the completeness of a cap ^ in 
PG(r, 2^). The key idea is making the more effective as possible the opera¬ 
tions that are performed more often. 

We identihed the following operations: 

1 . computing the points belonging to a line passing through two points of 

2. normalizing a point P hnding a representation with the leftmost non¬ 
zero coordinate equal to 0; 

3. keeping track if a point is saturated or not. 

1. is obtained representing the points of PG{r,2^) as integers, so that 
the sum of two points is computed as a bitwise exclusive or between integers. 
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Multiplication by a scalar for the points of ^ is precomputed at the beginning 
of the algorithm. 

2. is avoided by the way 3. is performed; we use a vector V of booleans 
that represents each point P of PG{r, 2^) 2^ — 1 times: we represent all vectors 
aP, a G F2fc\{0}; P is saturated is at least one of the aP is saturated. 

We used this algorithm for proving that the 95256-cap in PG{12, 4) of [13] 
is complete. 

When using this algorithm with greater values of r or k, it can became 
too memory consuming. In this case one either can reintroduce point 2. and 
represent in point 3. only the normalized version of P or can split V in 
subvectors Vi and test the covering of the points of Vi separately. In the 
latter way memory is saved, but computational time increases. However 
these computations can be performed in parallel. 
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